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QUANTUM  COMPUTING:  FROM  BRAGG  REFLECTIONS  TO  DECOHERENCE 

ESTIMATES 


Peter  Pfeifer  and  Chen  Hou 

Department  of  Physics,  University  of  Missouri 

Columbia,  MO  65211,  U.S.A. 

ABSTRACT 

We  give  an  exposition  of  the  principles  of  quantum  computing  (logic  gates,  exponential 
parallelism  from  polynomial  hardware,  fast  quantum  algorithms,  quantum  error  correction, 
hardware  requirements,  and  experimental  milestones).  A compact  description  of  the  quantum 
Fourier  transform  to  find  the  period  of  a function — the  key  step  in  Shot’s  factoring  algorithm — 
illustrates  how  parallel  state  evolution  along  many  classical  computational  paths  produces  fast 
algorithms  by  constructive  interference  similar  to  Bragg  reflections  in  x-ray  crystallography.  On 
the  hardware  side,  we  present  a new  method  to  estimate  critical  time  scales  for  the  operation  of  a 
quantum  computer.  We  derive  a universal  upper  bound  on  the  probability  of  a computation  to 
fail  due  to  decoherence  (entanglement  of  the  computer  with  the  environment),  as  a function  of 
time.  The  bound  is  parameter-free,  requiring  only  the  interaction  between  the  computer  and  the 
environment,  and  the  time-evolving  state  in  the  absence  of  any  interaction.  For  a simple  model 
we  find  that  the  bound  performs  well  and  decoherence  is  small  when  the  energy  of  the  computer 
state  is  large  compared  to  the  interaction  energy.  This  supports  a recent  estimate  of  minimum 
energy  requirements  for  quantum  computation. 


INTRODUCTION 

A quantum  computer  puts  to  use  the  fact  that  a quantum  particle,  such  as  an  electron  passing 
through  a double-slit  apparatus  or  an  atom  or  molecule  traversing  a matter-wave  interferometer 
[1],  can  be  in  two  different  locations  at  the  same  time.  By  equating  different  locations — as  a 
paradigm,  we  take  an  electron  in  the  lowest  orbit  or  in  an  excited  orbit  of  an  atom — to  binary 
digits  0 and  1,  one  may  interpret  the  time-evolving  state  of  the  particles  as  executing  several 
computations  at  the  same  time.  One  set  of  locations  at  a given  time  describes  the  result  of  one 
computation.  Thus  one  atom  can  do  two  computations  at  once;  two  atoms  can  do  four;  three 
atoms  can  do  eight.  The  challenge  is  to  coerce  the  atoms  to  follow  trajectories  that  amount  to 
meaningful  computations  and  to  read  out  a definite  result  from  the  multitude  of  computations 
occurring  in  parallel.  The  control  of  trajectories  is  the  hardware  part  of  the  challenge;  the  design 
of  useful  trajectories — algorithms  that  are  superior  to  classical  algorithms — is  the  software  part 
of  the  challenge. 

This  paper  is  a comprehensive  survey  of  both  aspects.  It  illustrates  in  terms  of  explicit  case 
studies  the  key  components  of  a quantum  computer,  namely,  the  design  of  a computational 
trajectory  (Shor’s  algorithm),  the  trajectory’s  superiorperformance  (exponential  parallelism), 
and  the  reliable  operation  of  the  computer  vis-a-vis  decoherencc,  with  emphasis  on  analogies  to 
familiar  physical  phenomena.  The  case  studies  apply  regardless  of  whether  the  computer  is 
implemented  in  terms  of  nuclear  spins,  ion  traps,  cavity  quantum  electrodynamics,  supercon- 
ducting currents,  or  spintronics.  General  references  are  [2-9]. 
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QUANTUM  VERSUS  CLASSICAL  COMPUTATION 


A classical  computer  manipulates  strings  of  N classical  bits,  (n,, with  -0  or  1 
(y»l.  in  such  a way  that  intermediate  states  of  the  computation  are  also  strings  of 
classical  bits.  A quantum  computer  manipulates  states  of  N two-level  atoms  or  any  other  two- 

level  entities,  |n, with  n.  «0  if  the  y-th  atom  is  in  the  ground  state  and  “1  if  it  is  in 

the  excited  state,  in  such  a way  that  intermediate  states  are  superpositions  of  the  states 
|rti,  The  2*'  states  |n,,  (“computational  basis”)  are  product  states  in  which  each 

atom  is  either  in  the  ground  state  or  excited  state,  and  rij  is  called  the  value  of  the  y'-th  qubit 
(“quanmm  bit");  they  represent  the  strings  of  classical  bits.  The  superpositions  include  states  in 
which  an  atom  no  longer  has  a sharp  value  of  (indefinite  bit  value),  and  states  in  which  an 
atom  no  longer  exists  in  a state  separate  from  the  other  atoms  (entangled  state);  both  have  no 
classical  counterpart. 

A quantum  computation  starts  with  a product  state  |n,,  lets  the  state  evolve  according 

to  the  Schrbdinger  equation,  r'h^|^(r))  - W(/)|tp(r)),  with  initial  condition  |tK0))  “ |" "w)  ^nd 

time-dependent  Hamiltonian  H{t)  driving  the  coupled  atoms,  and  ends  with  the  measurement  of 
the  values  of  the  qubits  of  the  state  |tp(r)).  The  Hamiltonian  generates  the  unitary  time-evolution 
operator  U (r)  which  takes  the  initial  state  into  the  final  state. 


IVKO)  - Tt%^-(ilfi)jH{s)dsy(Q))  - (/(r)|tp(0)) 


(1) 


with  T the  time-ordering  operator.  The  measurement  transforms  |V<(r))  into  the  output  state 

|n,' nj,)  with  probability  |(«',  ...,nJ^|V>(0)f-  The  output  is  probabilistic  because  quantum 

measurements  are  so.  The  computation  \n fi^)t-4|tp(t))  is  a unitary  transformation,  hence 

reversible;  the  readout  )^(r))i-»|n|' nj,)  is  a projection  (“collapse  of  the  wave  function"), 

hence  irreversible.  Thus,  to  perform  a specific  computation,  one  must  drive  the  atoms  with  a 
specific  Hamiltonian;  to  read  out  the  result,  one  mast  send  the  atoms  through  a series  of  state 
detectors. 

A quantum  computer  is  more  powerful  than  a cla.ssical  computer  for  two  reasons,  (i)  The 
quantum  state  space  is  much  larger  than  the  classical  state  space:  N qubits  can  be  in  an  infinite 
number  of  different  states  (any  point  on  the  unit  sphere  of  the  complex  Hilbert  space  spanned  by 
the  2*  basis  vectors  |«,,  N classical  bits  can  only  be  in  2*  different  states  (the  points 

where  the  2*  coordinate  axes  intersect  the  unit  sphere.  Fig.  1).  If  the  expansion  coefficients  of 
the  superpositions  can  be  controlled  with  accuracy  t (f  « 1),  the  sphere  Irosts  (Tfc''^""'’)  distinct 
states.  Thus  a quantum  computer  can  store  and  access  an  exponentially  large  number  of  states 
compared  to  a classical  computer,  (ii)  The  quantum  computer  operates  in  a massively  parallel 
way:  if  the  initial  state  is  the  uniform  superposition  of  all  basis  states, 

jV(0))  = 2-^'=  2 I", «,).  (2) 

» 

the  time  evolution  computes  simultaneously  t/(/)|n| for  all  2''  possible  inputs  |n n„) 

by  linearity  of  U{i)-  Hk  matrix  element  (n,' n;,|(/(r)|n,,  ...,«„)  is  the  probability  amplitude 

that  the  computation  converts  the  input  |n, n^}  into  the  output  |n,' rty),  along  all  possible 

classical  computational  paths  in  parallel  (Feynman's  path  integral).  A classical  computation  can 
follow  only  a single  path.  The  aim  is  to  choose  U{t)  so  that  computational  paths  of  no  interest 
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Figure  1.  A classical  bit  can  be  in  state  |0)  or  |1)  (left).  A quantum  bit  can  be  in  any  superposi- 
tion o|0)  + b|l),  where  the  complex  numbers  a and  b satisfy  (a|^  + |b|’  = 1 (right).  The  states  are 
shown  as  expectation  values  of  the  spin  vector. 

cancel  each  other  by  destructive  interference,  and  paths  of  interest  add  constructively  (selection 
of  relevant  computations  by  quantum  interference.  Fig.  2). 

Any  computation  f/(r),  also  called  a quantum  circuit,  can  be  approximated  by  sequential 
application  of  a finite  set  of  unitary  transformations  that  operate  on  only  one  or  two  qubits.  An 
example  for  such  a “universal  set  of  quantum  logic  gates”  is 

“(|0/X0y|H0r)(l/|H0(0/|-KK'/l)/^  (“Hadamardgate”),  (3) 

r-|0,.)(0,.|  + c‘”'^|l,.)(l,.|  CTgate”),  (4) 

-|0AX0Ah|0yl.X0;,|  + |l,0,)(l,4|  + |l,l4(l,^  (“controlled-notgate”),  (5) 

where  acts  only  on  qubit y,  and  In^n,  only  on  qubits y and  it.  They  correspond 

to  the  logic  gates  in  a classical  computer,  but  are  reversible  (the  classical  ‘and’  and  ‘exclusive  or’ 
gates  are  irreversible).  The  Hadamard  gate  transforms  the  states  |0^ ) and  |1^ ) into  the  superpo- 
sitions (|0j)±|l^))/-y2;  the  T gate  shifts  the  phase  of  the  excited  state  relative  to  the  ground 
state  by  ;r/4;  and  the  controlled-not  gate  flips  the  “target”  qubit  k if  and  only  if  the  “control” 
qubit y is  in  the  excited  state.  The  three  gates  are  the  analog  of  an  optical  beam  splitter,  phase 
shifter  (refractive  medium),  and  conditional  mirror,  respectively.  Only  the  Hadamard  gate  cre- 
ates multiple  computational  paths;  the  other  two  transform  a single  basis  state  into  a single  basis 
state.  To  approximate  a general  f/(r)  to  within  accuracy  e requires  OiA" N^[\n{A'* /e)]°) 
gates,  where  a « 2 (Solovay-Kitaev  theorem  [6])  and  the  leading  factor  4"  is  set  by  the  number 
of  matrix  elements  of  U(t).  For  special  computations,  however,  often  a much  smaller  number  of 
gates,  independent  of  e,  suffices,  as  we  shall  see  in  the  next  section. 


Figure  2.  Application  of  the  Hadamard  gate  to  state  |n  J creates  two  computational  paths,  one 
leading  to  jO^  ) with  probability  amplitude  1/^,  the  otner  to  |1^)  with  amplitude  (-1)'7a/2- 
Application  of  the  Hadamard  gate  twice  to  |0^ ) creates  four  paths,  leading  to  |0^ ) with  ampli- 
tude [(+lX+l)  + (+lX+ 1)1/2.  and  to  |1^)  with  amplitude  [(+!)(+ l)  + (+lX-l)1/2  (cancellation  of 
trajectories  leading  to  |I^  j). 


AN  EXPLICIT  EXAMPLE— FINDING  THE  PERIOD  OF  A FUNCTION 

A widely  used  class  of  encryption  systems  for  public  transmission  of  sensitive  data  (RSA 
public  key  cryptosystem  [10])  derive  their  security  from  the  difficulty  of  factoring  a large,  pub- 
licly transmitted  integer.  The  fastest  known  classical  algorithm  factors  an  N-bit  number  in  lime 
0(2“"'*''''’""^’’'').  Shot’s  celebrated  quantum  algorithm  [II,  12]  for  the  same  task  requires  only 
time  CXiV^jlnNXInlnN)).  The  key  task  in  Shot’s  algorithm  is  to  find  the  period  r of  a periodic 
function  /(jt)  related  to  the  number  to  be  factored,  and  the  tool  used  to  determine  the  period  is 
the  discrete  Fourier  transform.  The  values  of  the  Fourier  variable  at  which  the  transform  does 
not  vanish  yield  the  period,  just  as  in  x-ray  diffraaion  from  a crystal  the  scattered  wave  vectors 
at  which  the  x-ray  intensity  is  nonzero  yield  the  crystal  periodicity  (Bragg  reflection  at  wave 
vectors  equal  to  a reciprocal  lanice  vector). 

The  discrete  Fourier  transform  assigns  to  a given  complex-valued  function  fix)  (not  neces- 
sarily periodic)  the  function 


(6) 

where  both  x and y run  through  the  integers  0, 1 , . . 2"  - 1 

. Its  quantum  computation  reads 

>'v)"i/(f)2!LoVwki ■'^) ' 

(7) 

uup >v)> 

(8) 

x-x,2'*-'+jt,2*-^  + ..,  + jc^2". 

(9a) 

(9b) 

Equation  (8)  defines  the  operator  (/(r)  in  terms  of  the  computational  basis;  (9)  is  the  binary 
expansion  of  the  summation  variables  in  (7,  8);  and  Eq.  (7)  is  established  by  substitution  of  (6) 
and  (8).  The  sums  in  (7),  running  over  the  whole  basis,  completely  describe  the  function  and  its 
Fourier  transform.  The  implementation  of  (Jit)  in  terms  of  logic  gates  is  given  by 


I 


= 2''"^J][''_i(|0^)  + |i^ 

= /?(//, 5, A---S,^)(//,fl^3B„.--5,J---(//,,_A  -,x,).  (10) 


The  first  line  recalls  the  action  of  i/(t)  on  a basis  state;  the  second  line  shows  that  the  result  is  a 
product  state,  not  an  entangled  state  (the  binary  fractions  arise  from  the 

expansion  of  xy/2^);  and  the  third  line  expresses  the  product  in  terms  of  successive  applications 
of  the  Hadamard  gate,  the  controUed-phase-shift  gate. 


= Ko,)(0,0, 1+ 10^.1,)(0,1, 1+  |l,0,)(l,0* I + |l  .l,XlA 


(11) 


and  the  qubit  reversal  operator,  /?|n,,  (implemented  by  reading  the  qubits  in 

reverse  order),  on  the  original  basis  state. 

Equation  ( 10)  is  the  quantum  fast  Fourier  transform.  It  is  an  explicit  example  for  the  decom- 
position of  a quantum  circuit  into  a product  of  one-  and  two-qubit  operations  and  shows  that  the 
transform,  without  any  approximation,  can  be  carried  out  with  only  N{N  + 1)/2  operations.  This 
is  much  less  than  the  Solovay-Kitaev  bound  and  much  less  chan  the  0(2"  Al)  operations  in  the 
classical  fast  Fourier  transform.  The  N Hadamard  gates  aeate  2"  classical  paths  along  which 
the  state  |x,,  evolves  in  parallel. 

The  period  r is  obtained  by  measuring  the  qubits  of  state  (7).  The  measurement  yields  the 
output  state  |y,,  ....yj,)  with  probability 

|(v ,w|V'(0)f  = |/(y)f.  (12a) 

If  r divides  2",  the  probability  is 

l/(y)f-^|S;'./“'''’7(4-0  ify/2"-0,l/r,2/r (r-D/t,  (12b) 

|/(y)f-0  else,  (12c) 


from  (6).  Thus  as  advertised,  only  “good”  y’s,  related  to  r,  remain.  Outcomes  unrelated  to  r 
occur  with  probability  zero,  by  virtue  of  destructive  interference.  The  good  y’s  yield  r as  the 

largest  denominator  of  the  reduced  fractions  y/2",  as  y runs  through  all  outcomes  |y, y^). 

The  largest  denominator  from  m distinct  outcomes  (m  « r)  equals  r with  probability  better  than 
1 - (1-0.561 /Inin  r)"  (prime  number  theorem),  so  even  a period  as  large  as  2“"  is  found  with  a 
success  rate  better  than  99%  from  only  33  distinct  outcomes.  If  r does  not  divide  2",  destructive 
interference  is  not  perfect,  and  the  measurements  with  high  probability  yield  y’s  satisfying 
y/2"  -0, 1/r, 2/r, ...,  (r-l)/r  approximately.  In  this  case,  r is  found  from  a slightly  more 
involved  analysis,  but  still  with  only  0(lnN)  Fourier  transforms  and  measurements.  This  gives 
the  period  in  time  0{N^\nN).  With  this  efficiency,  it  turns  out  that  the  bottleneck  in  Shor’s 

algorithm  is  the  computation  of  f(x),  i.e.,  the  preparation  of  the  initial  state  2,  /(x)|x x,/), 

which  we  have  suppressed  here.  The  determination  of  the  period,  including  the  preparation  of 
the  initial  state,  requires  time  0(A'^(bN)(lnlnAl)). 
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Table  I.  Exponential  parallelism  of  the  quantum  Fourier  transform  from  polynomial  hardware. 
The  second  line  is  often  expressed  as  that  the  quantum  computer  can  do  2"  classical  computa- 
tions simultaneously.  The  abbreviations  in  the  last  column  stand  for  ‘effort’  and  ‘payoff. 


Number 

Logic  gates  acting  on  |V>(0)) 

A/(/V  + l)/2 

(10) 

Polynomial  e. 

Oassical  bit  strings  proces,sed  at  once  in  Ivlt)) 

2" 

(7) 

Exponential  p. 

Qassical  paths  followed  at  once  in  |(/r(t)) 

(10) 

Exponenbal  p. 

Distinct  constructive  interference  events  from  1^(0) 

2"r 

(12) 

Exponential  p. 

Repeated  computations  of  \y){t))  needed  to  find  r 

0(lnN) 

- 

Logarithmic  e. 

Table  I summarizes  the  algorithm  in  terms  of  parallel  computing  and  interference  of  compu- 
tational paths.  The  number  of  distinct  constructive  interference  events  increases  with  increasing 
period  r,  just  as  in  x-ray  diffraction  the  number  of  Bragg  reflections  increases  with  increasing 
size  of  the  unit  cell.  The  degree  of  entanglement  of  computational  states  varies  in  the  opposite 
direction:  If  / is  zero  everywhere  except  at  position  jt,  then  /has  period  r - 2*',  the  initial  state  is 
a product  state  by  |vt(0))  - 1 jt,,  ....  jr^ ),  and  so  is  the  final  state  (tp(/))  by  (10).  In  this  case,  neither 
the  initial  nw  final  state  is  entangled.  If  f(x)  - 2'"'^  for  all  x,  then  /has  period  r - 1,  the  initial 
state  equals  the  state  (2).  and  the  final  state  is  |vt(0)-|0,  ...,0).  In  this  case,  the  initial  state  is 
maximally  entangled,  and  the  final  state  is  not  entangled.  1^is  suggests  that  the  smaller  the 
period,  the  more  entangled  are  relevant  states  in  the  computation. 


OTHER  FAST  QUANTUM  ALGORITHMS 

Factoring  integers  can  be  reduced  to  finding  integer  solutions  of  Pell’s  equation,  x*  - cy  ^ - 1 , 
where  a is  a positive  non-square  integer,  but  a converse  is  not  known.  So  it  is  of  interest  that  a 
quantum  algorithm  exists  also  for  Pell’s  equation.  The  algorithm  is  due  to  Hallgren  [1 3]  and 
finds  the  solution  in  polynomial  time  as  in  Shot’s  algorithm,  instead  of  exponential  time  as  on  a 
classical  computer. 

Other  quantum  algorithms  that  outperform  classical  algorithms  by  orders  of  magnitude  are: 
Grover’s  algorithm  for  “finding  a needle  in  a haystack”  (search  of  an  item  in  a daubase  of  2" 
items  in  0(2'"^)  instead  of  0(2")  steps);  estimation  of  the  median  and  mean  of  2"  items  to 
precision  e in  0(N  /c)  instead  of  0{N  le^)  steps;  search  of  the  minimum  of  a function  sampled  at 
2"  points  in  0(2"'^)  instead  of  0(2")  steps;  search  of  two  distinct  pre-images  giving  the  same 
image  of  a two-to-one  function  sampled  at  2"  points,  in  0(2"'’)  instead  of  0(2"'^)  steps;  the 
Deutsch-Jozsa  algorithm  to  determine  if  2"  numbers  are  either  all  0 (“constant  function”),  or 
half  are  0 and  half  are  I (“balanced  function”),  in  one  instead  of  up  to  2"''  + 1 steps;  and  various 
allocation  tasks  and  game-theoretic  strategies. 

Quantum  algorithms  for  solving  complex  physical  problems  in  polynomial  time  instead  of 
exponential  time  include  [9]:  path  integration  with  respect  to  Wiener  measure;  quantum  random 
walks;  simulation  of  ±e  quantum  baker’s  map;  quantum  lattice-gas  model  for  the  time- 
dependent  Schrodinger  equation  of  many-body  systems;  and  the  dynamical  sign  problem  in 
fermionic  systems. 
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QUANTUM  ERROR  CORRECTION 


Noise  from  imperfect  computer  operation  poses  no  fundamental  barrier  to  large-scale  compu- 
tations. A quantum  error-correction  code  encodes  the  N “logical  qubits”  into  N'  “carrier  qubits” 
{N'  > N),  runs  the  carrier  qubits  through  a group  of  accordingly  encoded  logic  gates,  transforms 
the  noisy  state  by  appropriate  projection  operators  (“syndrome  measurements”)  and  unitary 
operators  (“recovery  of  original  carrier  qubits”)  into  an  error-corrected  state,  and  feeds  the  state 
into  the  next  group  of  gates.  Such  periodic  error  correction  prevents  accumulation  of  errors  in 
the  state  (“quantum  Zeno  effect”).  At  the  end  of  the  computation,  the  carrier  qubits  are  decoded. 
The  encoding  spreads  the  state  of  the  N logical  qubits  over  all  AT  carrier  qubits  so  that  when  the 
syndromes  are  measured,  no  information  about  the  state  of  the  logical  qubits  is  revealed  (“noise- 
less” or  “decoherence-free”  subspaces):  the  projections  preserve  superpositions  of  the  logical 
qubits,  and  the  state  of  the  carrier  qubits  is  highly  entangled  even  if  the  logical-qubit  state  is  not. 
Remarkably,  a discrete  set  of  corrections  can  correct  a continuum  of  errors.  For  a code  to  correct 
any  error  on  any  M carrier  qubits,  a necessary  condition  is 

N'fc4Af  + N (13) 

(Knill-Laflamme  bound),  and  a sufficient  condition  for  large  N is 

N IN' <1-2  [-xlog,  jt  - (I-  xllog^d  - Jt)]^.2„  (14) 

(Gilbert- Varshamov  bound).  E.g.,  a code  exists  which  encodes  one  logical  qubit  {N  - 1)  into  5 
carrier  qubits  (N'  - 5)  and  corrects  any  error  on  any  one  carrier  qubit  (M  1)  with  16  pairs  of 
syndrome  measurements  and  recovery  operations  (equality  in  the  IGiill-Laflamme  bound). 

Different  codes  require  different  encoding  of  gates,  and  of  interest  are  encodings  for  which 
an  error  in  an  input  carrier  qubit  or  the  gate  operation  propagates  only  to  a small  number  of 
output  carrier  qubits.  Specifically,  an  encoded  gate  is  called  fault-tolerant  if  a failure  with 
probability  p of  any  single  component  (e.g.,  one  of  the  10  “wires”  feeding  two  logical  qubits, 
each  encoded  by  5 carrier  qubits,  into  a controlled-not  gate)  introduces  an  error  in  two  or  more 
carrier  qubits  in  any  logical  output  qubit  with  probability  cp^  at  most,  with  c a constant  and  p 
small.  Such  a gate,  when  followed  by  error  correction  wifli  Af  - 1,  yields  an  error-free  output 
with  probability  1 - cp^  i.e.,  reduces  the  error  probability  from  p to  q>^  if  p<Mc.  Fault-tolerant 
Hadamard,  T,  and  controlled-not  gates  exist  with  c -10''  and  d«10^,  where  d is  the  number  of 
operations  on  carrier  qubits  needed  to  encode  and  error-correct  the  gate.  By  hierarchical  fault- 
tolerant  encoding  of  all  gates,  a computation  involving  L gates  can  be  carried  out  to  within 
accuracy  e using  only  0(Z.[ln(L/e)]'°‘’'^)  operations  on  carrier  qubits,  if  p<Mc  (threshold 
theorem  for  quanhim  computation).  Thus,  if  the  noise  in  individual  carrier  qubits  is  low  enough, 
p<l/c  -lO"^,  arbitrarily  large  computations  can  be  performed  because  the  overhead  for  error 
correction  grows  only  polynomial-logarithmically  with  the  size  of  die  computation,  L. 


EXPERIMENTAL  STATE  OF  THE  ART  AND  HARDWARE  REQUIREMENTS 

A remarkable  array  of  experimental  realizations  of  quantum  computing  devices  and  algo- 
rithms have  been  achieved  to  date.  Current  record  holders  with  respect  to  the  number  of  qubits 

i 

I 
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that  can  be  controlled  and  prepared  in  well-defined  states  include,  in  the  different  categories  of 
experimental  implementation: 

(i)  7-qubit  nuclear  magnetic  resonance  devices  (NMR);  three  'H  and  four  ”C  nuclei  in  trans- 
crotonic  acid,  each  in  the  spin  up  or  down  state,  with  a total  of  twelve  2-qubil  gates  driven  by 
radio-frequency  pulses,  preparing  the  “Schrtxlinger-cat  state”  (jOOOOOOO)  + 11 1 H 1 1 1))  and 
converting  it  into  the  state  |0000000)  [14];  and  five  '’F  and  two  ‘’C  nuclei  in  a custom-made 
molecule,  fully  implementing,  with  a sequence  of  300  radio-frequency  pulses,  Shor’s  algorithm 
to  factor  the  number  15  [15]; 

(ii)  a 4-qubit  ion-trap  device  (IT):  four  *Be*  ions  in  a linear  electromagnetic  trap,  each  in  one 
of  two  hype^e  Zeeman  levels  and  driven  by  Raman  transitions  [16]; 

(iii)  a 3-qubit  device  based  on  cavity  quantum  electrodynamics  (CQED):  three  Rb  atoms, 
each  in  one  of  two  Rydberg  states  (principal  quantum  numbers  49,  50,  and  51)  and  coupled  to  a 
cavity  mode  with  zero  or  one  photon,  driven  by  microwave  pulses  [17];  and 

(iv)  1 -qubit  devices  based  on  macroscopic  persistent-current  states  in  superconductors  [18], 

Computations  carried  out,  in  addition  to  Shot's  algorithm,  include:  Grover’s  algorithm  on  a 

2-qubit  NMR  device  (1998/99);  dynamics  of  quantum  harmonic  and  anharmonic  oscillators  on  a 
2-qubit  NMR  device  (1999);  the  Deutsch-Jozsa  algorithm  on  a 5-qubit  NMR  device  (2(XX3); 
finding  the  order  of  a permutation  on  a 5-qubit  NMR  device  (2000);  correction  of  any  one-qubit 
error  on  a 5-qubit  NMR  device  (2001);  noiseless  encoding  of  one  logical  qubit  in  a 3-qubit  NMR 
device  (2001);  quantum  lattice-gas  treatment  of  the  diffusion  equation  and  Burgers  equation  on  a 
2-qubit  NMR  device  (2002);  quantum  baker's  map  on  a 3-qubit  NMR  device  (2002);  protection 
of  an  IT  qubit  from  decoherence  (two  *Be*  ions  encoding  one  qubit  (2001));  controlled-not  gate 
on  a 2-qubit  IT  device  (one  *Be*  ion  implementing  a motional  and  an  internal-state  qubit  (2002)); 
and  Grover’s  algorithm  in  a 5-“qubit”  classical  optical  cavity,  demonstrating  that  the  algorithm 
requires  no  entanglement  (2(X)2). 

Any  experimental  realization  faces  three  challenges; 

(i)  Tlie  device  must  be  able  to  control  the  state  of  each  qubit  separately,  while  allowing 
neighboring  qubits  to  interact  with  each  other  for  the  operation  of  two-qubit  gates.  Control  of 
individual  qubits  is  achieved  by  different  chemical  shifts  in  NMR,  laser  beams  driving  ions 
spatially  separated  by  tens  of  pm  in  IT,  and  pulses  addressing  successive  atoms  traveling  at 
spatial  separation  of  several  cm  in  CQED.  Qubits  interact  via  spin-spin  coupling  in  NMR,  a 
shared  phonon  state  of  the  ions  (excitation  of  the  cenier-of-mass  miode)  in  IT,  and  electric  dipole 
coupling  between  each  atom  and  the  cavity  mode  in  CQED. 

(ii)  The  system  must  be  switchable,  so  that  interactions  executing  a prescribed  sequence  of 
gate  operations  (Hamiltonian  //(f))  can  be  turned  on  and  off  by  external  control.  Switching  is 
done  by  magnetic  field  pulses  in  NMR,  and  laser  pulses  in  IT  and  CQED. 

(iii)  The  computer  must  be  well  isolated  from  the  environment  so  that  the  decoherence  time 
r,,,  the  time  at  which  the  computer  and  the  environment  depart  significantly  from  a product  state 
(entanglement  of  the  computer  with  the  environment),  is  long  compared  to  f^,  the  time  it  takes  to 
operate  a single  gate.  At  time  f,,,  a generic  qubit  state  o{0)  + fi|l)  will  have  degraded  into  the 
mixture  |aP|0X0|-r|/>|^|l)(l|  (density  matrix),  which  no  longer  contains  the  interference  terms 
necessary  for  quantum  computation.  Good  isolation  is  provided  by  long  spin-spin  and  spin- 
lattice  relaxation  times  in  NMR  (f^  - 10"^  - 10*  s,  - lO"*  - 10"’  s),  long-liv^  hyperfine  levels 
and  stable  trap  and  laser  operation  in  IT  (r^  - 10''  - 10"  s,  » 10"’  - 10"’  s),  and  the  low  sponta- 
neous emission  rate  of  Rydberg  states  and  low  photon  escape  rate  from  the  cavity  in  CQED 
(f,-10-’-10"s,  f,-10-’-10^  s). 
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DECOHERENCE  ESTIMATES 


The  condition  is  a challenge  because  the  computer  must  be  a weakly  open  system  for 
time  /g,  but  a strongly  closed  system  for  time  t^.  The  condition  p < 1 /<•  in  the  threshold  theorem 
for  quantum  computation  is  even  more  demanding.  To  see  this,  let  |r/>j(/))  be  the  state  of  the 
computer  and  environment  interacting  with  each  other,  starting  from  a product  state  at  r = 0;  and 
let  |tl'„(0)  be  the  same  state,  but  noninteracting  (product  state  at  all  times); 


|Vi(0>  - rexp(-07fi)  + ff.(s)  + V(s)}dsjxi,M}VM), 

\ip„{t))  - |r exp|-(i7h)J^//,(s)*|]t//,(0))  |r exp|-(/7fi)J^//,(s)*j|Ve(0))j. 


(15) 

(16) 


Here  and  V'(r)  is  the  Hamiltonian  of  the  computer,  environment,  and  interaction 

between  the  two,  respectively;  and  |Vj(0))  and  is  the  initial  state  of  the  computer  and 

environment,  respectively.  Suppose  now  the  overlap  (“fidelity”), 

(17) 


decays  exponentially  with  time.  The  decay  constant  is  the  inverse  of  the  decoherence  time, 
f (f)  - , and  the  threshold  condition  requires  the  failure  probability  p to  satisfy 

p = <l/c,  (18) 


which  gives  the  condition  /(2c)  for  large  c.  Thus  the  computer  must  remain  isolated  from 
the  environment  for  about  lO*  gate  operations  for  sustained  computations.  Present  devices  can 
execute  about  300  operations  [IS]  and  are  far  from  this  goal. 

One  may  seek  to  reduce  and  use  (18)  to  estimate  how  short  should  be  for  the  threshold 
theorem  to  kick  in.  But  for  short  times,  the  fidelity  does  not  decay  exponentially,  and  Eq.  (18)  is 
not  appropriate.  Also,  the  fidelity  evolves  differently  if  a gate  is  switched  rapidly  instead  of 
slowly,  which  further  limits  the  validity  of  (18).  It  is  Aerefore  of  interest  that  a universal  bound 
for  the  fidelity  exists  which  suffers  from  none  of  these  limitations.  The  bound  reads 


F(/)  2 F.(f)  :=  cos |min|h"‘ A(s)<£r,  Jr/ljj, 

(19) 

Ms) - {%(s)y(s^^As)f , 

(20) 

valid  for  all  r.  It  requires  only  the  knowledge  of  the  product  state  (16)  and  interaction  V'(s').  The 
energy  uncertainty  A(i)  is  easy  to  evaluate  compared  to  the  task  of  computing  the  fidelity  from 
(15-17).  The  bound  is  a “worst-case”  decoherence  estimate:  if  FJj)  is  close  to  1,  so  is  F{t)  by 
f (f)  s 1 , and  decoherence  is  guaranteed  small.  F (r)  equals  1 at  f - 0,  decreases  with  increasing 
t (the  more  slowly,  the  weaker  the  interaction),  and  reaches  zero  (trivial  bound)  when  the  integral 
exceeds  jdi/2.  The  time  t at  which  F_(r)  has  dropped  to  the  value  1/e,  fi''/JA(j)</j-=1.19, 
gives  a lower  bound  for  the  decoheience  time,  t s (Fig.  3).  A particularly  appealing  form  of 
inequality  (19)  is  the  upper  bound 
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p{t)  s sin^|min|ft'‘ &(s)ds, 


for  the  probability  of  the  quantum  computation  |V'c(0)  (first  factor  in  (16))  to  “fail"  after  time  r 
due  to  imperfect  isolation,  p(r) 1 - F^(t).  The  term  “fail"  is  short  for  the  interacting  state  to 
give  possibly  different  results  than  the  noninteracting  state.  The  bound  (21)  estimates  the  failure 
probability  rigorously,  at  any  point  in  time,  and  without  reference  to  gate-switching  and  decoher- 
ence times  (parameter-free  decoherence  estimate),  unlike  (18). 

The  result  (19)  is  an  application  of  the  variational  principle  of  Pfeifer  and  Frohlich  [19]  (see 
also  [20])  which  states  that  if  |v'(0)  is  a solution  of  the  time -dependent  Schrbdinger  equation 
with  Hamiltonian  H(t)  and  |v,(t))  is  a trial  state  with  |tp,(0))-|tp(0)),  then 

Kv»(t)|v;,(0)|acos|min|h''j^  |l-|V',(r))(t/',(s)|)(W(j)-ih^)|t//,(s))|di,  ar/2|j.  (22) 

The  norm  under  the  integral  tracks  the  amount  by  which  the  trial  state  does  not  satisfy  the  Schro- 
dinger  equation.  For  the  choice  //(r)-//^(/)  + //.(O  + Vf/),  |t/K0)”|V’,(0),  and  |V't(0}”|V'n(0). 
the  norm  reduces  to  A(r),  whence  (19). 

To  explore  the  performance  of  (19),  we  have  computed  (22)  numerically  for  a model  system. 
The  system  is  a particle  in  a step  potential  in  one  dimension,  which  splits  an  incoming  wave 
packet  into  a transmitted  and  reflected  wave  packet,  similar  to  a Hadamard  gate.  But  here  the 
potential  models  the  interaction  of  the  computer  and  environment.  The  initial  state  is  a Gaussian 
wave  packet.  The  trial  state  is  a Gaussian  wave  packet  with  position  and  momentum  equal  to  the 
classical  motion  across  the  potential  step,  neglects  the  reflected  wave  packet,  but  includes  the 
spreading  of  the  transmitted  packet  (“quasi-free  panicle”).  Table  II  maps  the  computer  dynamics 
onto  the  wave-packet  dynamics.  Figure  4 shows  the  result  for  incident  kinetic  energy  100  times 
the  energy  of  the  potential  step,  and  initial  wave-packet  width  1/10  of  the  initial  distance  from 
the  step.  It  is  a sample  from  15  computations,  varying  in  initial  energy  and  width.  When  the 
wave  packet  reaches  the  step  and  splits  into  a large-amplitude  transmitted  and  small-amplitude 
reflected  component,  the  overlap  drops  from  1.0(K)  to  0.997;  when  the  transmitted  component  is 
well  separated  from  the  reflected  component,  the  overlap  rises  back  to  1.000.  In  comparison,  the 


Figure  3.  Lower  bound  t for  the  decoherence  time  (schematic).  The  upper  curve  is  FU). 
assumed  to  obey  F(t)-e''"‘  for  large  t (definition  of  r^);  the  lower  curve  is  f (/). 
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Table  II.  Wave-packet  model  (WP)  of  the  dynamics  of  a quantum  computer  (QC). 


oc 

WP 

Isolated  computer; 

• Hamiltonian 

• State 

Wc(0  + Wc(0 

\wM 

Quasi-firee  particle 
|V.(0) 

Nonisolated  computer; 

• Hamiltonian 

• State 

H,(.t)  + H,(.t)  + V(t) 

im) 

Particle  in  step  potential 

m 

Fidelity; 

• Lower  bound 

• Lower  bound  = 1 if 

• Lower  bound  - 1 if 

Eq.  (19) 

V(f)  = 0 

Energy  of  l^nlO))  is  large 

Eq.(22) 

Potential  step  = 0 
Traversal  of  step  is  fast 

lower  bound  drops  to  0.95  when  the  wave  packet  reaches  the  step,  and  remains  at  that  value 
because  the  integral  in  (22)  accumulates  all  departures  of  the  trial  state  from  the  exact  state;  a 
temporary  departure  produces  a permanent  drop.  The  bound  does  not  drop  further  because  the 
trial  state  is  asymptotically  exact  on  the  far  side  of  the  step.  Thus  the  bound  performs  very  well. 
It  performs  well  whenever  the  wave  packet  traverses  the  step  fast,  i.e.,  the  arrival  time  at  the  step 
is  long  and  the  residence  time  at  the  step  is  short. 

Fast  traversal  of  the  step  translates  into  a large  energy  of  the  computer  state  compared  to  the 
interaction  energy.  Thus  the  probability  of  a gate  to  fail  at  time  f^,  p(/j),  is  small  if  the  energy 
supplied  to  switch  the  gate,  E^,  is  large;  (dpIdE^),  <0.  This  agrees  with  the  recent  estimate, 
p - h/(t^E^),  of  the  smallest  achievable  failure  probability,  or  minimum  energy  required  to 
switch  the  gate,  under  isolation  [21].  Our  result  is  a fuat  step  to  extend  the  estimate  to  the  case 
where  the  computer  is  not  isolated. 


Figure  4.  Numerical  evaluation  of  inequality  (22).  Left;  overlap  of  the  exact  state  and  the  trial 
state  [left-hand  side  of  (22)].  Right;  lower  bound  of  the  overlap  [right-hand  side  of  (22)].  Time 
is  in  units  of  h/(energy  of  the  potential  step). 
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CONCLUSIONS 


We  have  reviewed  quantum  computing  as  an  interplay  of  deterministic  time  evolution  and 
probabilistic  measurement  outcomes,  constructive  and  destructive  interference,  open  and  closed 
system  dynamics,  and  theory  and  experiment.  We  have  developed  an  inequality,  Eq.  (21).  which 
decomposes  the  failure  probability  of  the  computer,  due  to  imperfect  isolation,  into  a part  that 
maintains  coherence,  |V'„(r)),  and  the  part  that  causes  dccohcrence,  V(s).  The  inequality  reduces 
the  estimation  of  failure  probabilities  to  a computationally  easy  problem  for  a wide  range  of 
compniter  designs.  In  an  exploratory  model  computation,  we  find  that  the  failure  probability  can 
be  made  small  by  making  the  energy  stored  in  the  computer  large.  This  suggests  there  may  be  a 
variety  of  ways,  to  be  uncovered  yet,  to  minimize  decoherence  at  given  V (s). 
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